Interaction Effects in the Mesoscopic Regime: 
A Quantum Monte Carlo Study of Irregular Quantum Dots 
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We address the issue of accurately treating interaction effects in the mesoscopic regime by investigating the 
ground state properties of isolated irregular quantum dots. Quantum Monte Carlo techniques are used to calcu- 
late the distributions of ground state spin and addition energy. We find a reduced probability of high spin and a 
somewhat larger even/odd alternation in the addition energy from quantum Monte Carlo than in local spin den- 
sity functional theory. In both approaches, the even/odd effect gets smaller with increasing number of electrons, 
contrary to the theoretical understanding of large dots. We argue that the local spin density approximation over 
predicts the effects of interactions in quantum dots. 
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PACS numbers: 73.23.Hk, 73.63.Kv, 02.70.Ss 

The interplay between correlations and quantum mechani- 
cal interference of electron states, long at center stage in con- 
densed matter physics, has been traditionally investigated in 
disordered systems but can also be probed in confined sys- 
tems, such as quantum dots. 1,2 In these latter, the confinement 
leads to mesoscopic fluctuations 1 - 3 which in turn modify the 
role of Coulomb repulsion between electrons within the dots. 
Quantum dots offer the great practical advantage of experi- 
mental tunability in the study of this interplay*^ 

Quantum dots of different size give rise naturally to differ- 
ent descriptions. For small dots, including both vertical! and 
few electron lateral 5 dots, circular symmetry is preserved and 
plays a critical role. In this limit quantum Monte Carlo (QMC) 
calculations have been performed4i2i2ii£ as well as simpler 
density functional theory (DFT) simulations. 1112 A compar- 
ative study for the weak interaction regime 7 ' 8 confirmed the 
validity of the DFT method in this small dot limit. 

For large irregular dots, 1,3 on the other hand, all spatial 
symmetries are broken. For a sufficiently irregular shape, the 
motion of electrons within the dot is chaotic, which then jus- 
tifies modeling the single-particle energies and wave func- 
tions by random matrix theory and random plane waves, 
respectivelyj^ii Furthermore, interaction effects in these 
larger dots are often treated within the random phase ap- 
proximation (RPA) for gas parameter r s (the ratio between 
the interaction energy and kinetic energy, formally defined as 
r s = 1/ y / 7mao for 2D bulk systems, which thus identifies 
the strength of the interaction) of order 1 or smaller i 2 * 1 ^ The 
"universal Hamiltonian" picture 2 - 16 that emerges leads to sta- 
tistical predictions for various quantities, such as the ground 
state spin or addition energy. 1718 Once temperature is taken 
into account, 18 these are in good agreement with experimental 
data^-* 2 ^-^-* 2 - One notable feature is a substantial difference at 
zero temperature between dots containing an odd number of 
electrons, N, and those in which this number is even. Experi- 
ments have not to date performed at a sufficiently low temper- 
ature to probe this feature. This even/odd effect persists for 
large N with an essentially unchanged magnitude provided 
that r s remains constant. 

In order to go beyond statistical predictions and address 



features of individual irregular dots, an approach which ac- 
curately treats the combination of mesoscopic fluctuations 
and interaction effects is needed. DFT appeared a natural 
choice for such studies, and, indeed, microscopic calculations 
of ground state energies for large irregular dots (N ~ 200) 
were carried out within the framework of the local spin den- 
sity approximation (LSDA)^ 3 * 2 !^ The statistics of the LSDA 
results turned out, however, to be in qualitative disagreement 
with the earlier predictions, even for the modest interaction 
strengths (r s ~ 1.5) that are experimentally relevant; for in- 
stance, in the LSDA results at zero temperature, the even/odd 
effect is nearly absent. In fact, there were several indications 
of stronger interactions in LSDA than those obtained from 
RPA. The striking discrepancy between the two approaches - 
both of which are believed to be valid in the range of r s con- 
sidered - combined with the absence of experimental statistics 
for low temperature keeps this problem open. 

Here we take up the issue of accurately treating interac- 
tion effects in the mesoscopic regime. We consider irregular 
quantum dots with up to 30 electrons. The lack of symmetry- 
induced shell structure makes irregular dots qualitatively dif- 
ferent from circular dots; in particular, the mesoscopic inter- 
ference effects are both more subtle and more generic. In this 
regime, where the universal Hamiltonian picture is not ex- 
pected to hold because of the modest size, we use QMC to 
treat the interactions much more carefully than in LSDA. To 
this end, we present QMC calculations of the addition energy 
and ground state spin for such dots, and compare to corre- 
sponding LSDA results. 

We consider a model quantum dot consisting of electrons 
moving in a two dimensional plane, with kinetic energy 
(— i V?)> and interacting with each other by long-range 
Coulomb repulsion Q3i<i \ r i ~~ r il _1 )- All energies are ex- 



= 1, with 



pressed in atomic units, defined by h = e /e = 
electronic charge e, effective mass m*, and dielectric constant 
of the medium e. The electrons are confined by an external 
(quartic) potential 



V ext (x,y) = a 
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where r — y/x 2 + y 2 . This simple form of V cx t breaks all 
symmetries except time reversal invariance. It leads to chaotic 
motion of the electrons inside the dot, 26 which is the experi- 
mental situation for large dotSpiA2i 2 £i 2 ii2 2 i in fact, not only the 
bare V CK t but also the self-consistent potential leads to chaotic 
dynamicsi2i2i We have studied potential Q for a range of pa- 
rameters and report here results for a — 0.002 (which controls 
r s ), b = 7r/4, 7 between 0.1 and 0.2 (break spatial symme- 
tries), and A between 0.53 and 0.67. For these parameters, 
the dynamics in the bare potential is chaotic. We accumulate 
statistics for six dots formed by different sets of parameters (A 
and 7) from the above range. We study dots with N = 10 to 
30 electrons which yields a range r s — 1.8 to 1.3i2I 

Variational (VMC) and diffusion (DMC) Monte Carlo 
techniques 28 were used to calculate the energies E(N, S) of 
our model quantum dots for each N and spin S. We investi- 
gated S = 0, 1, and 2 for even N, and S = 1 /2, 3/2, and 5/2 
for odd N. For a given V ex u me ground state energy Eqs and 
the ground state spin Sgs were determined for each N. 

The trial wave function used in QMC, ^t> is written as 
a linear combination of products of up- and down-spin Slater 
determinants multiplied by a Jastrow factor. Each Slater deter- 
minant is constructed from single-particle Kohn-Sham (KS) 
orbitals obtained using the LSDA functional. The Jastrow 
factor effectively describes the dynamic correlation between 
the electrons coming from their mutual repulsion, whereas 
the near-degeneracy or static correlation is taken into account 
by having more than one determinant. We optimize the Jas- 
trow parameters and determinant coefficients by minimizing 
the variance of the local energy 



In a second stage, we use fixed-node DM 1 



2S.30 



to project 



the optimized many-body wavefunction onto a better approx- 
imation of the true ground state. The fixed-node DMC en- 
ergy is an upper bound to the true energy and depends only 
on the nodes of the trial wave function, i.e., only on the linear 
combination of determinants. (The Jastrow factor affects the 
statistical error of the energy but not its expectation value.) 
The statistical error in the energy E'dmc (N, S) obtained in 
this way is smaller than the single-particle mean level spacing 
A by about two orders of magnitude and hence is insignifi- 
cant. The systematic error from the fixed-node approximation 
of the many-body wave function is, however, difficult to esti- 
mate though experience suggests it is often small. We have 
included Slater determinants for which the sum of the KS 
single-particle energies are up to A greater than the sum of 
the KS single-particle energies for the ground state KS de- 
terminant. This amounts to taking mostly one and sometimes 
two or three Slater determinants in the *S?t expansion. Increas- 
ing the energy window from A to 3A failed to reduce -Edmc 
although it sometimes reduced TSvmc- 

We present the distribution of Sgs obtained from both 
DMC and LSDA -in Fig. 1(a). Within a model of effectively 
non-interacting electrons, Sgs is or 1/2 for even TV and 
odd N, respectively, due to standard up/down filling of the or- 
bitals. We see that the probability of finding a non-trivial Sgs, 
i.e. not zero or half, is substantial. Interestingly, this probabil- 
ity is reduced in DMC calculations compared to LSDA. Note 
that the differences between the two distributions, although 
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FIG. 1 : (a) Distribution of Sgs from DMC and LSDA calculations. 
Shaded histograms are for dots with even N while unfilled bars are 
for odd N; data is collected for N = 10-30 and six realizations of 
Vixt, which are also used to estimate the statistical error. Though 
the differences are small, LSDA predicts a larger probability of non- 
trivial Sgs- (b) Distribution of the difference in "spin gap" obtained 
using DMC and LSDA, normalized by the mean level spacing A. 
The large width of the distribution (crp) indicates a significant differ- 
ence between the two techniques. Note that x (defined in the figure) 
is primarily positive; negative values occur predominantly when the 
ground state has non-trivial spin. LSDA is, therefore, making non- 
trivial spin states more probable by lowering their energies compared 
to DMC results. [A sliding Gaussian window of width 0.08 is used 
to give a smooth estimate for P(x).] 



clearly visible, are not large and are therefore not much big- 
ger than the statistical error given the relatively small data set 
(21 x 6 = 126 cases total). There are, however, significant cor- 
relations between the LSDA and DMC results. In fact, Sqs 
from DMC is, up to one exception, always the same as or 
lower than that from LSDA. As a consequence the statistical 
error on (S^f DA - S^f ) is only 30% of its value. Thus, 
there is a clear difference between LSDA and DMC in the 
predicted ground state spin. 

The ground state spin shows the difference between LSDA 
and DMC results at only a coarse level. To obtain a more 
detailed understanding, we focus on the "spin gap", 5, which 
we define as 5 = E(S = 1) - E(S = 0) for even N and 5 = 
E(S = 3/2) - E(S = 1/2) for odd N. Thus S is the amount by 
which the higher spin state differs in energy from the lower. 
Changing the spin of a dot from or 1/2 to a higher value 
involves a competition between the single-particle energy cost 
and the exchange energy, —JS(S + 1), gainiS 

The key result for comparing DMC and LSDA is shown in 
Fig. 1(b): the distribution of (<5dmc — <5lsda)/A. Note, first, 
that the distribution is broad (full width ~ 0.4A); it is of order 
the energy required to flip Sgs from to 1 for an even dot, 
A — 2 J, assuming a realistic value of the exchange parameter 
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FIG. 2: Top: The addition energy from both DMC and LSDA for 
one realization of V cx t as a function of the number of electrons on 
the dot. Note the large mesoscopic fluctuations in A 2 E around the 
expected overall decrease. Bottom: Fluctuation in the addition en- 
ergy after removing the smooth part, normalized to the mean level 
spacing A. We see that the fluctuations are of order A and that they 
are somewhat larger for DMC results than for LSDA. 




FIG. 3: Distribution of the normalized fluctuations in the conduc- 
tance peak spacing, s = (A 2 E - (A 2 E)tit)/A from DMC (top) and 
LSDA (bottom) calculations. The DMC distributions for even (solid) 
and odd (dashed) N are quite different while there is less difference 
for LSDA. The standard deviation in the different cases, a, quantifies 
this contrast. (A sliding Gaussian window of width 0.17 (for even) 
and 0. 1 (for odd) is used to obtain a smooth curve.) 



J~ 0.35 for r s ~ 1.5. Second, note that the spin gap in DMC 
tends to be larger than that in LSDA. This indicates that the 
strength of interactions in LSDA is overestimated. Finally, we 
have studied this quantity separately for the smaller (N= 10- 
20) and larger (N — 20-30) dots. We have not found any size 
dependence - results in both ranges of N are the same as in 
Fig. 1(b) within our statistical accuracy. These observations, 
together with the results of Fig. 1(a), show that LSDA unduly 
favors non-trivial spin states. 

The ground state spin distribution has implications for the 
distribution of the spacing between Coulomb blockade con- 
ductance peaks, 17 through its relation to Eqs- In the nearly 
isolated dot limit, the spacing between the Coulomb blockade 
conductance peaks is proportional to the the addition energy 1 
denned by 

A 2 E(N) = E GS {N + 1) + E GS (N-1)-2E GS (N) . (2) 



For noninteracting electrons one would have 
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for even N 
for odd N 



(3) 



where e< are the energies of the single-particle states. Note 
the sharply different characteristics of even and odd N dots. 
Interactions reduce this strong even/odd effect. 

The behaviorof A 2 E(N) for a particular realization of V cx t 
is presented in the top panel of Fig. 2. Similar qualitative 
behavior is observed for all other configurations we studied. 
An overall decrease of A 2 E with N is expected due to the 
increase of effective capacitance of the dot and, hence, the 



decrease in the classical charging energy as the dot gets big- 
ger. On top of the mean behavior, we clearly see strong meso- 
scopic fluctuations arising from the interplay of electron inter- 
action and interference effects in the irregular dots. The fluc- 
tuations seem to be slightly larger in the DMC results than in 
LSDA. To focus on these fluctuations, we subtract the smooth 
classical part using a linear fit, and present the fluctuating part 
normalized by the mean level spacing A (which is the natural 
scale of these mesoscopic fluctuations) in the lower panel. 

To get a more quantitative picture, we plot the distribution 
of the normalized addition energy fluctuations in Fig. 3. Note 
the larger difference between the two distributions obtained 
from DMC than between those obtained from LSDA. In the 
DMC results, the width of the distribution for even N is sig- 
nificantly larger than that for odd dots, and the even N distri- 
bution has a long tail reminiscent of the Wigner surmise for 
the distribution of £jv/2 _e w/2-i found using random matrix 
theory (RMT) [see Eq. 

Looking at the data for smaller and larger dots separately, 
we find a strong trend shown in Fig. 4: the even/odd effect 
in the DMC data decreases significantly as N increases. If 
extrapolated to much larger N, this trend, also present in 
our LSDA results, contradicts the predictioniil& of combining 
RPA interactions with an RMT treatment of the single particle 
statistics. 

In conclusion, we have used quantum Monte Carlo to ac- 
curately investigate the role of interactions in the mesoscopic 
regime. We find that for irregular dots with a gas parameter 
r s ~ 1.5 and electron number in the range 10 to 30, DMC 
calculations show (1) mesoscopic fluctuations of the addition 
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energy, (2) a substantial probability of non-trivial ground state 
spin, and (3) a significant even/odd effect in the addition en- 
ergies. In comparison to LSDA, DMC typically predicts a 
somewhat larger spin gap; as a consequence, it has a tendency 
to find smaller ground state spins and a somewhat stronger 
even/odd effect in the addition spectra. These findings sug- 
gest that LSDA, as compared to DMC, in some sense over 
predicts the effect of interactions. It is interesting to note that 
a similar conclusion concerning overly strong interactions in 
LSDA was reached in the large dot regime 25 using the Struti- 
nsky analysis scheme. 



FIG. 4: Distribution of the normalized fluctuations in the conduc- 
tance peak spacing, from DMC calculations in the range of N — 10- 
20 (left) and N = 20-30 (right). For smaller N, the even/odd effect 
is stronger: cr cvcn = 0.46 ± 0.08, while cr od d =0.18 ± 0.02. On the 
other hand for larger N, the even/odd effect is significantly reduced: 
o-ovcn = 0.30 ± 0.04 and a od d = 0.25 ± 0.04. A similar qualita- 
tive trend is found in the LSDA results, though quantitatively it is 
somewhat weaker. 
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